%% --------- SCRIPT EJEMPLO NEWTON-RAPHSON PARA TP1 ------------------

close all  %cierra todas las ventanas
clear all  %elimina variables del entorno

a = [3;1;0];
b = [1;1;0];

x = [2;3;0];                     %ejemplo con 1 piedra
% x = [ [2;3;0] , [2;0.5;0] ];   %ejemplo con 2 piedras

t = 0 : 0.05 : 1;

C = A(0.5,a,b,x);

plot( t, C - A(t,a,b,x), '.-' )

x0 = 0.1; tol = 0.000001; maxnumit = 100;

xfinal = newton_raphson(x0,a,b,x,C,tol,maxnumit)

%lo mismo con precision simple:
% xfinal = newton_raphson(single(x0),single(a),single(b),...
%     single(x),single(C),single(tol),maxnumit);